SUPERNOVA REMNANTS IN THE SEDOV EXPANSION 
PHASE: THERMAL X-RAY EMISSION 
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ABSTRACT 

Improved calculations of X-ray spectra for supernova remnants (SNRs) in the Sedov- Taylor 
phase are reported, which for the first time include reliable atomic data for Fe L-shell lines. This 
new set of Sedov models also allows for a partial collisionless heating of electrons at the blast wave 
and for energy transfer from ions to electrons through Coulomb collisions. X-ray emission calcu- 
lations are based on the updated Hamilton-Sarazin spectral model. The calculated X-ray spectra 
are succesfully interpreted in terms of three distribution functions: the electron temperature 
and ionization timescale distributions, and the ionization timescale averaged electron tempera- 
ture distribution. The comparison of Sedov models with a frequently used single noncquilibrium 
ionization (NEI) timescale model reveals that this simple model is generally not an appropriate 
approximation to X-ray spectra of SNRs. We find instead that plane-parallel shocks provide a 
useful approximation to X-ray spectra of SNRs, particularly for young SNRs. Sedov X-ray models 
described here, together with simpler plane shock and single ionization timescale models, have 
been implemented as standard models in the widely used XSPEC vll spectral software package. 



Subject headings: ISM: supernova remnants - X-rays: ISM 



INTRODUCTION 



X-ray astronomy has advanced significantly in 
recent years because of improved X-ray instrumen- 
tation. Several recent X-ray satellites provided 
complementary capabilities: high spatial resolu- 
tion and mapping of the entire sky (ROSAT), 
spatially-resolved spectroscopy (ASCA), and broad- 
band coverage (BeppoSAX). In the field of su- 
pernova remnants (SNRs), this has led to many 
advances, such as the discovery of the nonther- 
mal synchrotron X-ray emission produced by en- 
ergetic electrons accelerated at a shock front, 
identification of a number of new neutron stars 
and associated synchrotron nebulae, discovery of 
new SNRs, and identification of young, ejecta- 
dominated SNRs. Spatially- resolved X-ray spec- 
troscopy, provided by the ASCA satellite, played 
a crucial role in these discoveries. The ASCA 
database now contains a large number of excel- 
lent spectra of SNRs. Observations with the next 
generation of X-ray satellites ( Chandra and XMM- 



Newton) have begun to provide even better quality 
data, ensuring continuing progress in the SNR re- 
search. 

While the new X-ray observations led to signif- 
icant progress, X-ray spectral data on SNRs, and 
ASCA data in particular, have been mostly under- 
utilized. The fundamental reason for this unsat- 
isfactory situation is that X-ray emitting plasmas 
in SNRs are not in ionization equilibrium, signif- 
icantly complicating analysis of their X-ray spec- 
tra. Until recently most X-ray observers have not 
had access even to simple nonequilibrium ioniza- 
tion (NEI) models. The simplest NEI model con- 
sists of an impulsively heated, uniform and homo- 
geneous gas, initially cold and neutral. While this 
constant temperature, single-ionization timescale 
model is often used in analysis of X-ray spec- 
tra of SNRs, in most cases it is unlikely to be 
a good approximation to shock-heated plasmas 
in SNRs. Because X-ray spectra are sensitive to 
the detailed structure of a SNR, the best descrip- 
tion of SNR spectra is through hydrodynamical 
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modeling, coupled with X-ray emission calcula- 
tions. This approach is obviously impractical in 
most cases. That leaves idealized hydrodynami- 
cal structures, such as self-similar hydrodynami- 
cal solutions, as the best choice for general NEI 
models. The Sedov- Taylor self-similar solution, 
describing a point blast explosion in a uniform am- 
bient medium, is the most well-known and useful 
of such structures. 

X-ray emission calculations for the Sedov mod- 
els have been performed numerous times in the 
past, with the most extensive set of calculations 
presented by Hamilton, Sarazin, & Chevalier 
(1983). The shape of the X-ray spectrum, spa- 
tially integrated over the whole remnant, depends 
on the shock speed v s , the ionization timescale 
r (which we define as the product of the rem- 
nant's age to and the postshock electron number 
density n es ), and the extent of electron heating at 
the shock front. In addition, as typical for most 
X-ray emitting thermal plasmas, line strengths 
and to a lesser extent continua depend on heavy- 
element abundances. Sedov models have been re- 
calculated using updated atomic data (Kaastra & 
Jansen 1993) and successfully applied to modeling 
of X-ray spectra of SNRs in the Large Magellanic 
Cloud (Hughes, Hayashi, & Koyama 1998). 

The most uncertain aspect of X-ray modeling 
in the framework of the Sedov model is the extent 
of electron heating in a blast wave. Because of the 
low densities typical for the interstellar medium 
(ISM), shocks in the ISM are collisionlcss, where 
the ordered ion kinetic energy is dissipated into 
random thermal motions through collisionless in- 
teractions involving magnetic fields. Electrons, 
whose collisions with ions produce the observed 
X-ray emission from SNRs, are most likely also 
heated to some degree in collisionlcss shocks, but 
the details of this process are not well under- 
stood. In the extreme case of a very high Mach- 
number shock, in the shock frame both electrons 
and ions enter with the shock velocity v s , so that 
after randomization one expects kTi ~ rriiVg and 
kT e ~ m e v 2 s , that is, T e /Ti ~ 1/1836, unless some 
collisionless heating of electrons (at the expense of 
ions) also occurs. Hamilton et al. (1983) consid- 
ered two extreme possibilities, either full electron- 
ion equilibration at the shock front, or electron 
heating only through Coulomb collisions in the 
postshock gas, without any collisionless heating at 



the front itself. Hughes et al. (1998) used these 
two limiting classes of models in their modeling 
of LMC SNRs, and found better (but not per- 
fect) agreement for models without any collision- 
less heating. While these results suggest that the 
full electron-ion cquipartition might not hold in 
SNRs, the extent of collisionlcss electron heating 
remains an open question. 

These questions are critical for observations 
since both the X-ray continuum and lines origi- 
nate from electrons: from electron bremsstrahhmg 
or electron impact excitation of ions. Until X-ray 
observations are able to resolve thermal Dopplcr 
widths of emission lines, ion energy remains invis- 
ible except as a pool from which to heat electrons. 
The detailed shape of the electron distribution af- 
fects both continuum and line emission. This in- 
troduces another potential problem: the possible 
deviation of the electron energy distribution func- 
tion from a Maxwcllian distribution, in particular 
the presence of a high-energy tail. Electrons in 
such a tail might modify the high-energy X-ray 
spectrum of a SNR through interactions with ions 
and magnetic fields. 

Theoretical modeling of collisionlcss shocks, 
performed with the purpose of accounting for the 
observed intensity and properties of Galactic cos- 
mic rays, revealed that energetic particles gener- 
ated within the shock should strongly modify the 
shock structure. Such cosmic-ray modified shocks 
transfer a significant amount of kinetic energy into 
cosmic rays, which then diffuse upstream of the 
original shock and partially decelerate incoming 
gas in a region ahead of it (the shock precursor). 
The temperature of the thermal gas is much lower 
than in a standard (nonmodificd) shock with the 
same velocity, because of this efficient transfer of 
kinetic energy into cosmic rays. This will affect X- 
rays produced by thermal plasmas in cosmic-ray 
modified shocks, because of strong dependence of 
X-ray emission on gas temperature. The modeling 
of modified shocks has been usually done in the 
framework of steady-state, plane-parallel shocks, 
preventing a direct comparison with the Sedov 
model discussed above. Such a comparison is now 
possible, as time-dependent, cosmic-ray modified 
shock models have been constructed recently by 
Berezhko, Yelshin, & Kscnofontov (1994, 1996). 

The continuing improvement in theory just de- 
scribed is not followed by a similar improvement 
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in the X-ray data analysis. Instead of using Se- 
dov models (and even more sophisticated hydrody- 
namical models, including those with cosmic-ray 
modified shocks), many observers still rely on the 
ionization equilibrium models, which arc generally 
not acceptable for SNRs. Sometimes spectra are 
modeled by Gaussian lines and free-free continua, 
a purely phenomenological approach which does 
not allow for determination of physical parameters 
in X-ray emitting plasma. The frequently used 
constant temperature, single ionization timescale 
NEI model is better than an equilibrium plasma 
model, but it is still not appropriate for SNRs. 
In order to make full advantage of current and 
future X-ray data on SNRs. a concentrated ef- 
fort is required to close the gap between the rela- 
tively sophisticated theory and the current crude 
data analysis. The theoretical work described here 
focuses on a most basic set of models which is 
necessary for analysis of X-ray data, and which 
should be routinely used by X-ray observers. We 
present new calculations of spatially-integrated X- 
ray emission from Sedov models in § 2, allowing 
for partial electron heating at a shock front. A 
preliminary version of this work was presented in 
Lyerly et al. (1997). In § 3, wc discuss our results 
in terms of suitably chosen temperature and ion- 
ization timescale distribution functions, and com- 
pare Sedov models with simpler NEI models such 
as a single timescale, constant temperature model 
or a plane-parallel shock. The sensitivity of X-ray 
spectra to model parameters is discussed in § 4. 
We present a summary of our results in § 5. 

2. X-RAY EMISSION CALCULATIONS 
2.1. Sedov- Taylor Dynamics 

We base our spectral calculations on a Sedov- 
Taylor (ST) blast wave model. The ST model 
begins with a spherically symmetric point blast 
expanding adiabatically into a uniform ambient 
medium with the adiabatic index 7 = 5/3. Only 
the shocked ISM is considered in the ST model. 
Therefore, the shocked ISM mass must signifi- 
cantly outweigh the mass of the material ejected 
by the supernova (SN), usually at an age of a 
few thousand years. The governing parameters 
of this phase of SNR evolution are the initial en- 
ergy of the blast E, the density of the ambient 
medium p , and the time since the explosion t - 



These parameters can be easily arranged to form 
a unit of length, related to the shock radius by 

r s = 1.15 (^q/po) 1 ^ 5 - These three initial parame- 
ters completely determine the dynamics of a Sedov 
model. For example, we might choose an initial en- 
ergy of 10 51 ergs, an upstream hydrogen number 
density n = 1 cm -3 , and a time i of 10 11 sec- 
onds, yielding a radius of 2.5 x 10 19 cm. We now 
know the shock speed, v s — 0.4r s /i , and also the 
postshock values of pressure, density, and velocity 
through the Rankine-Hugoniot shock conditions. 
The self-similar analytical solution (Sedov 1946; 
Taylor 1950) gives the pressure, density and veloc- 
ity profiles behind the shock, and it allows us to 
find pressure and density variations as a function 
of time for each fluid element. This information is 
necessary for determination of the ionization state 
for each ion species. 

Hamilton et al. (1983) showed that the 
shape of the X-ray spectrum of a spatially- 
integrated Sedov model does not depend on all 
three parameters mentioned above, but only on 
shock velocity v s (or alternatively shock temper- 
ature T s = 1.37 x 10 7 (V1000 km s_1 ) 2 K = 
1.18(^/1000 km s -1 ) 2 keV) and parameter rj = 
n^E. (We have assumed cosmic-abundance 
plasma, where electrons are provided by fully ion- 
ized H and He, with a negligible contribution from 
heavier elements. The mean mass per particle \x 
is then equal to 0.6 in units of m p .) In order to 
simplify comparison with other NEI models, we 
use ionization timescale To instead of 7/, where 
we define tq as the product of the postshock 
electron density n es and the remnant's age to. 
Note that for cosmic-abundance plasma and the 
strong shock Rankine-Hugoniot jump conditions, 
n es = 4.8rto- From equation (4b) in Hamilton 
et al. (1983), ionization timescale To is equal to 
6.03 x lO 11 ^ 3 ^ 5/6 cm" 3 s, where 7751 = ry/lO 51 
ergs cm~ 6 and T 7 = T s /10 7 K. X-ray spectra de- 
pend also on electron temperature T e which is 
discussed next. 

2.2. Electron Temperature 

The mean gas temperature T in Sedov models is 
plotted in Figure 1 as a function of the normalized 
emission measure EM(r) = J^n e n p dr/ J^ s n e n p dr 
(Because we assumed normal abundances, with a 
cosmic H/He ratio and a negligible contribution 
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Fig. 1. — Electron temperature T e vs emission 
measure EM (both normalized to 1 at the shock) , 
for Sedov model with T s = 3 x 10 7 K, r = 5 x 10 11 
cm~ 3 s, and with full ion-electron equipartition 
(T e = T; solid curve), no collisionless heating of 
electron at the shock (dotted curve), and partial 
electron heating at the shock with the postshock 
electron temperature T es = 0.5T S = 1.5 x 10 7 K 
(dashed curve). 
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Fig. 2. — Ratio of electron temperature T e to 
the mean gas temperature T vs emission measure 
EM, for the Sedov model with the same T s and To 
as in Figure 1 , but with various amounts of colli- 
sionless heating ((3 ranges from 0.0 at the bottom 
to 0.9 at the top). 



to electron density from elements heavier than He, 
models and calculations discussed in this work are 
not applicable to heavy-element plasmas where 
these assumptions break down.) T is equal to the 
shock temperature T s at the shock (EM = 1), 
and increases toward the remnant's interior. But 
for calculations of X-ray spectra we need to know 
the electron temperature T e , which in general is 
not equal to either the mean temperature T or the 
ion temperature T^. These three temperature are 
equal only if electrons and ions are equilibrated at 
the shock front through plasma instabilities gen- 
erated within the collisionless plasma. Because 
kinetic energy is initially transferred to thermal 
motions of ions in collisionless shocks, full equili- 
bration demands the presence of a very efficient 
mechanism for transfer of energy from ions to 
electrons in collisionless shocks, at all shock Mach 
numbers. This is unlikely based on available ob- 
servational and theoretical evidence. For exam- 
ple, observations of ultaviolet (UV) lines from a 
fast nonradiative shock in SN 1006 (Laming et al. 
1996) suggest that there is little energy transfer 
from ions to electrons in shocks where neutral H 
is present ahead of the shock front. In an extreme 
case of a very inefficient heating of electrons at 
the shock, the postshock electron temperature T es 
is much less than the postshock ion temperature 
Ti S , and the ratio T es /Ti S (or T es /T s ) is close to 
0. Electrons are then heated mostly in Coulomb 
collisions with ions in the hot shocked plasma 
downstream of the shock front. The most general 
situation, most likely to happen in SNRs shocks, 
is that electrons are partially heated in collision- 
less shocks and additional Coulomb heating occurs 
downstream of the shock. 

The extent of electron heating in shocks can 
be parametrized by the temperature ratio (3 = 
T es /T s , which can vary from nearly in shocks 
with negligible collisionless heating to 1 when this 
heating is very efficient. In general, (3 might de- 
pend on the shock Mach number and other shock 
parameters, possibly decreasing with increasing 
shock Mach number (Bykov & Uvarov 1999). This 
inverse correlation with the shock Mach number 
is supported by most recent observational evi- 
dence based on optical observations of nonradia- 
tive shocks (Ghavamian 1999; Ghavamian et al. 
2000). Because of the poor understanding of elec- 
tron heating in collisionless shocks, which may de- 
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pcnd not only on the shock Mach number but also 
on other shock parameters such as the magnetic 
field orientation, we neglect dependence of on 
shock parameters, i. e., we assume that does 
not vary during the SNR evolution. The valid- 
ity of this assumption is questionable, and should 
be examined in more detail once good quality 
spatially-resolved X-ray spectra of SNRs become 
available. Once a value of is chosen, we can cal- 
culate the electron temperature T e in Sedov mod- 
els using the procedure described by Borkowski, 
Sarazin, & Blondin (1994), which takes into ac- 
count Coulomb heating in the shocked plasma. 
(The formula (D2) in this paper contains a ty- 
pographical error; the second term on the right 
hand side of this expression should be replaced by 
2.51n((l + /3 1 /2)/(i _ /3-1/2)) _ spi/Ztf + 3 )/ 3 .) 

The electron temperature T e of a fluid element de- 
pends on its location within the remnant, on T es 
and T s (hence on 0), and on the SNR ionization 
timescale r . 

The extreme case of inefficient electron heating 
at the shock is shown in Figure 1 for a Sedov model 
with T s = 3 x 10 7 K and T es = K (0 = 0), and 
To = 5 x 10 11 cm~ 3 s. The electron temperature 
rises from at the shock (EM = 1) to an asymp- 
totic finite value at the SNR center (EM = 0), the 
behaviour already noted by Itoh (1978, 1979) and 
Gronenschild & Mewe (1979) in their numerical 
simulations and discussed in more detail by Cox 
& Anderson (1982) and Hamilton et al. (1983). 
This is in contrast to very high electron tempera- 
tures in the interior of models with efficient heat- 
ing (T es = T s or = 1). It is obvious that mod- 
els with inefficient electron heating differ greatly 
from models with efficient heating. Another solu- 
tion, with intermediate — 0.5 (T es — 0.5T S ) and 
the same t = 5 x 10 11 cm -3 s, is also shown in 
Figure 1. 

The relative importance of Coulomb heating 
and collisionlcss heating in the Sedov model just 
discussed can be inferred from Figure 2, where 
we plot the ratio of electron temperature T e to 
the mean temperature T for a number of differ- 
ent values of parameter 0. For large 0, the curves 
are nearly straight parallel lines, with relatively 
unimportant Coulomb heating. The electron tem- 
perature profiles in these models arc nearly the 
same as in the Sedov model with 0=1, except 
for the overall temperature scaling which is pro- 



portional to 0. At low 0, the Coulomb heating 
dominates, resulting in highly curved lines. The 
efficiency of the Coulomb heating increases with 
increasing nAi/T 3 / 2 (Itoh 1978), where n is the 
total particle density of a fluid element, At is the 
time elapsed since this element was shocked by 
the blast wave, and T is the mean temperature of 
this fluid clement. At the blast wave and in the 
remnant's interior, nAt/T 3 ^ 2 is small, so that the 
Coulomb heating is less efficient there than inside 
the swept shell. This is why the curves in Figure 2 
attain a maximum at EM ~ 0.5 and drop to for 
EM = and 1. Note that because of the scaling 
just mentioned, these curves are identical for all 
Sedov models with the same value of tq/T^ 2 . 

The temperature profiles presented in Figures 1 
and 2 are merely examples, with a fixed value of 
tq/T^ 2 , which sample a small range of the pos- 
sible parameter space (0, r /T s 3 ^ 2 ). The effects 
of varying r /T s 3 ^ 2 are shown in Figure 3 for the 
same parameter r] as in Figures 1 and 2, but at dif- 
ferent times in the remnant's evolution (1/4, 1/2, 
1, 2, and 4 of 5 x 10 11 cm~ 3 s). The require- 
ment of constant r\ demands that T s is equal to 
3 x 10 7 K/(n t/5 x 10 11 cur 3 s) 6 / 5 . Two cases of 
= and 0.5 are shown. At early times, the fast 
collisionless heating in a model with = 0.5 dom- 
inates the Coulomb heating, leading to very differ- 
ent electron temperature profiles than for = 0. 
This is the case discussed until now. The temper- 
ature profiles become similar at later stages in the 
remnant's evolution, depending rather weakly on 
0, because of the increased role of the Coulomb 
energy transfer. At r = 20 x 10 11 cm~ 3 s (curves 
at the top of Figure 3), noticeable differences are 
seen only in the remnant's interior, where densi- 
ties are low and temperatures are high, leading to 
a long electron-ion equlibration time scale. Note 
again that these temperature profiles are identi- 
cal for all Sedov models with the same values of 
to/T 3/2 and 0. 

2.3. Ionization Calculations 

The ionization state of each fluid element de- 
pends on its ionization timescale t = J*° n e dt, 
where t s is the time when this fluid element was 
shocked, and on how the electron temperature T e 
varied as a function of its ionization time scale, 
from to r, during the remnant's evolution. The 
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Fig. 3. — Ratio of electron temperature T e to 
the mean gas temperature T vs emission measure 
EM, for the Scdov model with 77 = 8.88 x 10 51 
ergs cm" 6 , at five different times (increasing from 
bottom to top). Curves are labeled by the post- 
shock temperature. Two cases with (3 = (solid 
curves) and (3 — 0.5 (dotted curves) are shown. 



0.40 




0.10 7 



0.00 E , , , 1 

0.001 0.010 0.100 1.000 

EM 

Fig. 4. — Normalized ionization timcscale r vs 
normalized emission measure EM for Sedov mod- 
els. 



normalized ionization timescale t/to is plotted in 
Figure 4. It starts at zero at the shock (EM = 1), 
then increases linearly in the swept SNR shell, but 
in the remnant's interior it is again small because 
of low gas densities there. Note that the distri- 
bution of t is quite broad, and is poorly repre- 
sented by a S function (the approximation involved 
in using a single-ionization-timescale NEI model). 
The most highly ionized atomic species are usually 
found near the maximum of r, while the least ion- 
ized material is found just behind the shock and 
in the remnant's interior. We determined the elec- 
tron temperature as a function of r for each fluid 
element following the procedure described above. 

We tracked the ionic states in each fluid element 
by solving the time-dependent ionization equa- 
tions for each abundant element (C, N, O, Ne, Mg, 
Si, S, Ca, Fe, and Ni). We used an eigenfunction 
method to solve these equations efficiently (Masai 
1984, Hughes & Helfand 1985, Kaastra & Jansen 
1993, Borkowski et al. 1994). Physical processes 
included in these calculations are collisional ion- 
izations by electrons, both direct and involving 
excitations followed by autoionizations, and ra- 
diative and dielectronic recombinations. We used 
collisional ionization rates from Arnaud & Ray- 
mond (1992) for Fe, and from Lennon et al. (1987) 
for other elements. We also updated recombina- 
tion rates in the Hamilton & Sarazin (1984) X-ray 
code. Radiative recombination rates were taken 
from Pcquignot, Petitjean, & Boisson (1991), and 
rates for dielectronic recombination towards He- 
like, Li-like, and Be-like ions from Hahn (1993), 
and Tcng, Xu, & Zhang (1994a, b). Recombina- 
tion rates for Fe ions were taken from Arnaud & 
Raymond (1992). Except for H and He, all ele- 
ments are assumed to be neutral initially, a good 
approximation under most circumstances because 
of the insensitivity of X-ray spectra to the detailed 
ionization state of preshock plasma. 

2.4. Spectral Code 

We used an updated version of the X-ray code 
written by Hamilton & Sarazin (1984) to calcu- 
late X-ray spectra. We updated collision strengths 
for H-likc ions using fits provided by Itikawa et 
al. (1985) for 0+ 7 and by Callaway (1994) for 
other H-like ions. Collision strengths for He- 
like ions were taken from Sampson et al. (1983), 
Keenan et al. (1987), and Kato & Nakazaki (1989). 
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We also updated atomic data for fluorescent Fc 
Ka lines produced by inner-shell ionization. K- 
shell collisional ionization rates by electrons for 
all Fe ions were calculated according to Beigman, 
Shevclko, & Tawara (1996), and then combined 
with the fluorescence yields of Kaastra & Mewe 
(1993) to obtain strengths of fluorescent Ka lines 
for Fe +0 - Fe+ 21 . The fluorescent yield for Fc+ 22 
was taken from theoretical calculations of tran- 
sition probabilities and autoionization rates by 
Seely Feldman, & Safronova (1986). Energies, 
transition probabilities, and excitation rates by 
electrons for Fe L-shell transitions for Fe +16 - 

Fe +23 

were kindly provided to us by Duane Liedahl 
(see Liedahl, Osterheld, & Goldstein 1995 and 
Mewe, Kaastra, & Liedahl 1995 for the description 
of these theoretical calculations and their imple- 
mentation into the equilibrium ionization thermal 
MEKAL model). 

2.5. Sedov Models in XSPEC 

Modeling X-ray spectra of SNRs within the 
framework of Sedov models demands an efficient 
interface between observations and the calculated 
model spectra. Such an interface is provided by 
the well-known and widely-used software package 
XSPEC (Arnaud 1996). We wrote efficient FOR- 
TRAN programs for calculations of Sedov X-ray 
spectra, using methods and techniques just de- 
scribed, and constructed Sedov models in XSPEC 
v9 and vlO. (The most recent version of these 
models was implemented by Keith Arnaud into 
XSPEC vll.) The model parameters are: post- 
shock temperature T s , post-shock electron tem- 
perature T e , ionization timescale To, and elemen- 
tal abundances. Following standard practice used 
with thermal models in XSPEC, there are two ver- 
sions of Sedov models, where heavy-element abun- 
dances may be varied individually or together. 

Calculations of Sedov models are computation- 
ally intensive, resulting in very long computa- 
tional times when X-ray data sets are fitted with 
these models. An even more troublesome aspect is 
the lack of smooth convergence when fitting with 
Sedov models. Keith Arnaud solved this prob- 
lem by modifying a routine in the FTOOLS soft- 
ware package which was used to tabulate MEKAL 
model. With this routine we can now tabulate 
spatially-integrated Sedov models in a FITS file, 
which can then be used in conjunction with a stan- 



dard XSPEC model "atable" to make fits to the 
data. This procedure is necessary for the robust 
performance of the Sedov models in XSPEC. 

2.6. Calculated Spectra: Effects of Colli- 
sionless Electron Heating 

The flexible computer interface provided by the 
XSPEC package allows us to easily generate Sedov 
X-ray spectra in the whole parameter space rele- 
vant for SNRs, in the energy range from 0.1 keV to 
30 keV, and to simulate observations using spec- 
tral response matrices and effective areas of mod- 
ern X-ray satellites. There are clearly a number of 
interesting issues associated the use of Sedov mod- 
els, a lot of them discussed in detail by Hamilton 
ct al. (1983), which can be efficiently examined 
within the XSPEC package. Here we restrict our- 
selves to the effects of collisionless electron heating 
on X-ray spectra of SNRs. 

X-ray spectra for Sedov models with n = 8.88 x 
10 51 ergs cm~ 3 are shown in Figure 5 at differ- 
ent times during the remnant evolution, and at 
each time for three values of (3: 0.0, 0.5, and 1.0. 
All spectra were folded through the pre-launch in- 
strumental response of the back-illuminated ACIS 
CCD detector (S3) onboard the Chandra Obser- 
vatory, and include variations of the telescope's 
effective area with the photon energy. The second 
set of spectra from the top of Figure 5 corresponds 
to models whose electron temperature profiles are 
shown in Figure 1. These spectra differ substan- 
tially from each other, and become harder with the 
increasing amount of collisionless heating. Partic- 
ularly large differences are seen at high energies, 
reflecting large differences in temperature profiles 
seen in Figure 1. It should be obvious that, in 
general, it is necessary to consider not only the ex- 
tremes of very efficient and completely inefficient 
electron heating, but also models with partial elec- 
tron heating at the shock. 

Differences between X-ray spectra of Sedov 
models with different amounts of collisionless heat- 
ing become smaller as the remnant ages, because 
of the increased role of the Coulomb heating in 
older remnants without full electron-ion equilibra- 
tion at the shock front (see Figure 3). This is 
clearly seen in Figure 5, where there is hardly any 
difference between the coolest models with and 
without collisionless heating. Most of the X-ray 
emitting plasma in these models is actually not 
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very far from coronal equilibrium ionization. 
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Fig. 5. — X-ray spectra for Sedov models with 
T] = 8.88 x 10 51 ergs cm~ 6 and with various 
amounts of collisionlcss heating: full electron-ion 
cquipartition (solid curve), no collisionlcss heating 
of electrons at the shock ( dashed curve) , and par- 
tial electron heating at the shock with the post- 
shock electron temperature T es = 0.5T S (dotted 
curve) . Models are displaced from each other by a 
factor of 10 in the vertical direction and they are 
labeled by the postshock temperature (in keV). 



3. INTERPRETATION OF X-RAY SPEC- 
TRA 

X-ray and UV spectra of stellar coronae showed 
the presence of multi-temperature hot plasmas, 
which led researchers to consider temperature 
distribution functions in their efforts to under- 
stand the observed coronal spectra. This approach 
proved quite useful in analysis of X-ray spectra of 
stellar coronae. It would be advantageous to have 
a similar framework available for SNRs, so we gen- 
eralized the idea of the temperature distribution 
function used for plasmas in ionization equlibrium 
to NEI Sedov models. 

3.1. Temperature and Ionization-Timescale 
Distribution Functions 

The spatially-integrated X-ray spectra for Se- 
dov models can be succesfully interpreted in terms 
of the following three distributions: the electron 
temperature, ionization timescale, and ionization 
timescale-averaged electron temperature distribu- 
tion functions. We have already shown the ion- 
ization timescale distribution in Figure 4, and ex- 
amples of the electron temperature distribution in 
Figure 1. These two distributions do not contain 
enough information to determine ionic states, be- 
cause of the adiabatic expansion of the gas follow- 
ing its passage through the blast wave. This ex- 
pansion results in a decreasing mean temperature 
for each fluid element. This is demonstrated in 
Figure 6, where we plot temperature variations for 
several fluid elements in the Sedov model. Partic- 
ularly large variations are seen for fluid elements 
which were shocked early in the remnant's evo- 
lution. Another reason for electron temperature 
variations is energy transfer from ions to electrons 
in models with unequal ion and electron tempera- 
ture. In principle a detailed account of the electron 
temperature variations is needed for each fluid el- 
ement, just as we did for Sedov models, but a sim- 
pler approach in terms of an average temperature 
provides a better insight while maintaining good 
accuracy. A suitable average electron temperature 
is the ionization timescale-averaged temperature 
(T e ), defined as f*° T e (t)n e (t)dt/r. In Figure 7, 
we plot the ratio (T e )/T e for Sedov models with 
T s = 3 x 10 7 K, t = 5 x 10 11 cm" 3 s, and vari- 
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ous amounts of collisionlcss electron heating. For 
equal ion and electron temperatures this ratio is 
significantly greater than unity in the remnant's 
interior because of the adiabatic expansion. In 
the opposite limit of no collisionlcss heating, (T e ) 
is slightly less than T e in the postshock region, 
because the efficient transfer of energy from ions 
to electrons more than compensates for adiabatic 
expansion. However, even in this case the temper- 
ature ratio becomes larger than unity in the rem- 
nant's interior, demonstrating the effectiveness of 
the adiabatic expansion. After passing through 
the shock front, electrons in a fluid element are 
heated by Coulomb collisions with ions, but be- 
cause of the adiabatic expansion the efficiency of 
Coulomb heating decreases, and the electron tem- 
perature attains a maximum and then starts to 
decline with time. Even modest amounts of colli- 
sionless heating affect the ionization timescale av- 
eraged distribution function in the remnant's inte- 
rior and at the shock front (Fig. 7), where the ion- 
ization timescale is low and the Coulomb heating 
is least effective (Fig. 2). As the collisionless heat- 
ing becomes more important, the ratio (T e )/T e 
increases until it reaches its upper limit in fully 
equilibrated models. 

The temperature profiles presented in Figure 7 
are again merely examples, with a fixed value of 
tq/Ts^ 2 , which sample a small range of the pos- 
sible parameter space (/3, t /Ts^ 2 ). The effects 
of varying t /Ts^ 2 are shown in Figure 8, for Se- 
dov models with electron temperature distribu- 
tion functions shown in Figure 3. The ionization 
timescale averaged distributions differ from each 
other, reflecting the relative importance of adia- 
batic cooling and Coulomb heating at various lo- 
cations within the remnant. However, they are 
all bounded from below by models with low ef- 
ficiency of Coulomb heating and no collisionlcss 
heating, and from above by models with equal ion 
and electron temperatures. 

The two temperatures, the electron tempera- 
ture T e and the average electron temperature (T e ), 
may be used to obtain an approximate functional 
relationship of electron temperature with ioniza- 
tion timescale for each fluid element. A linear re- 
lationship is simple and sufficient, starting with 
temperature 2(T e ) — T e at the shock (ionization 
timescale equal to 0), and ending with T e at the 
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Fig. 6. — Temperature of fluid elements in the 
Sedov model, normalized to their final values, vs 
normalized ionization timescale t/tq. The tem- 
perature history is shown for 4 different fluid ele- 
ments (solid curves), labeled by their normalized 
emission measure EM. A linear approximation 
discussed in the text is also shown (dotted curves). 
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Fig. 7. — Ratio of ionization timescale-averaged 
electron temperature (T e ) to final electron T e vs 
emission measure EM for Sedov models shown in 
Figures 1 and 2, and with various amounts of col- 
lisionless heating. Curves are labeled by (3. 
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Fig. 8. — Ratio of ionization timescale-averaged 
electron temperature (T e ) to final electron T e vs 
emission measure EM for Scdov models shown in 
Figure 5, with r\ = 8.88 x 10 51 ergs. Two cases with 
(3 = (solid curves) and (3 — 0.5 (dotted curves) 
are shown. Curves are labeled by the shock tem- 
perature. Dashed curve corresponds to model with 
full electron-ion equilibration. 



current local ionization timescale r. We show this 
linear approximation in Figure 6 for a model with 
equal ion and electron temperatures. The agree- 
ment with the exact Sedov solution is excellent 
for the bulk of the X-ray emitting material, al- 
though it becomes somewhat less accurate in the 
remnant's interior where temperature variations 
are more pronounced. This approximate descrip- 
tion of electron temperature variations allows us to 
calculate ionization fractions and then X-ray spec- 
tra, and compare them with the exact X-ray cal- 
culations for Sedov models described in § 2.6. We 
find excellent agreement between the two sets of 
calculations for models with equal ion and electron 
temperatures, with typical differences at most of 
order of few percent (for cosmic abundance Sedov 
models at the ASCA or Chandra CCD spectral res- 
olution) . The agreement is slightly worse for mod- 
els without any collisionless heating, but the dif- 
ferences are rarely at the ~ 10% level and they are 
limited to a few line complexes at lower energies. 
Because such differences are smaller than current 
uncertainties in atomic data, we did not find it 
necessary to develop a more accurate approxima- 
tion to describe the evolution of electron temper- 
ature in a fluid element. The current linear ap- 
proximation is clearly sufficient for calculations of 
X-ray emission from the Sedov blastwave. While 
it offers negligible computational advantages over 
exact calculatons for Sedov models, we find it use- 
ful for interpretation of X-ray spectra. This ap- 
proximate description may also be very convenient 
in more complex situations, which might require 
multi-dimensional hydrodynamical simulations. 

X-ray spectra of SNRs are apparently described 
very well by the electron temperature, ionization 
timescale, and ionization timescale averaged elec- 
tron temperature distribution functions. High- 
quality X-ray data may allow us to determine 
these distribution functions without reliance on 
detailed hydrodynamical calculations. This is 
analogous to determinations of differential emis- 
sion measure in stellar coronae, without reliance 
on coronal models. The experience gained in anal- 
ysis of stellar coronae suggests that such an em- 
pirical approach should be useful in analysis of 
X-ray spectra of SNRs. But instead of one tem- 
perature distribution function needed in a stellar 
corona, three distribution functions are required 
to describe NEI X-ray spectra of SNRs. A purely 
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empirical determination of these distribution func- 
tions might then be a difficult task, so we expect 
that its combination with hydrodynamical model- 
ing will prove most useful in the future. 

3.2. Constant Temperature, Single Ioniza- 
tion Timescale Models 

The most common NEI model in use today 
is the constant temperature, single ionization 
timescale model. In terms of the distribution 
functions just discussed, electron temperature and 
average electron temperature are assumed con- 
stant, and the ionization distribution shown in 
Figure 4 is approximated by a delta function. As 
already noted by researchers (e. g., Masai 1994), 
the latter approximation is particularly trouble- 
some. If it is used to describe a plane-parallel 
shock or a Sedov model, prominent lines produced 
by strongly underionized plasma immediately be- 
hind the shock will be missing from the model 
spectrum. This is demonstrated in Figure 9, 
where we again show the calculated X-ray spec- 
tra for Sedov model with solar abundances, and 
with the same parameters as in Figures 1 and 7 
(T s = 3 x 10 7 K = 2.59 keV, t = 5x 10 11 cm" 3 
s). We fitted a simple NEI model just described 
to these spectra, using the XSPEC software pack- 
age, and the best fits are shown in Figure 9. The 
very poor match between the calculated spectra 
and the best fit model spectra demonstrates that 
the model under consideration is not adequate to 
represent X-ray spectra of SNRs. The current 
practice of using the constant temperature, single 
ionization timescale NEI model to represent X- 
ray spectra of SNRs is not valid in general. This 
simple NEI model may be used judiciously under 
special circumstances, e. g., for an isolated dense 
cloud of gas overrun by a shock long time ago or 
for SN ejecta in old remnants which were thermal- 
ized early in the remnant's evolution. However, 
these are exceptions to the rule because X-ray 
emission in SNRs is generally produced behind 
shocks, much as in Sedov models. 

3.3. Plane-Parallel Shock Models 

The failure of a constant temperature, single 
ionization timescale model is caused by an im- 
proper approximation of the ionization timescale 
distribution (Fig. 4) by the delta function. A 
broad distribution such as seen in Figure 4 is 
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PHOTON ENERGY (keV) 



Fig. 9. — Fits to X-ray spectra of Sedov mod- 
els of Figure 3, for clarity displaced from each 
other by factor of 10 and shown by solid curves 
(top - T es = T s , middle - T es = 0.5T S , bottom 
- T es =0). A constant temperature, single ion- 
ization timescale model (dotted curves) provides a 
poor fit when compared with plane-parallel shock 
(dashed curves). 
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needed in order to provide a good approximation 
to Sedov models. A simple and obvious choice is 
a plane-parallel shock characterized by a constant 
postshock electron temperature T e and by its ion- 
ization age t s — n e t s , where n e is the postshock 
electron density and t s is the shock age. This 
model differs from a single ionization timescalc 
NEI model in its linear distribution of ionization 
timescale vs emission measure (at the shock the 
ionization timescale is 0, and it attains the max- 
imum t s for the material shocked earliest). Such 
a simple shock provides a much better approxi- 
mation to the Sedov models shown in Figure 9, 
where we also plot the best fit shock model. The 
difference between Sedov models and the shock 
model is noticeable only at low energies, particu- 
larly for models with no collisionless electron heat- 
ing. For equal ion and electron temperatures, final 
values of the free parameters of the fit give a very 
good approximation to the mean electron temper- 
ature (T e — 3.26 keV vs the emission- weighted 
electron temperature of 1.27T S = 3.29 keV in the 
Sedov model) and the mean ionization timescale 
(r s /2 = 9.88 x 10 10 cm" 3 s vs (r) = 0.202r = 
1.01 x 10 11 cm~ 3 s, where (r) is the emission- 
measure averaged ionization timescale in Sedov 
models). Reasonable agreement is also obtained 
for the other two noncquipartition models with 
T es = 0.5T S = 1.30 keV and T es = 0: shock tem- 
peratures 2.15 keV and 1.58 keV vs the emission- 
weighted electron temperatures of 0.74T S = 1.91 
keV and 0.53T S = 1.37 keV in Sedov models, and 
t s /2 = 7.84 x 10 10 cm- 3 s, 6.50 x 10 10 cur 3 s 
vs 1.01 x 10 11 cm~ 3 s in Sedov models. This 
example demonstrates that plane-parallel shocks 
might serve as useful approximations to Sedov 
models with high electron temperatures in the 
Chandra spectral range. From the computational 
point of view, plane-parallel shocks with constant 
electron temperature require little computational 
overhead when compared with the single ioniza- 
tion timescale model (when ionic states are cal- 
culated with the eigenfunction technique). But 
unlike the latter model, plane-parallel shocks are 
useful for fitting X-ray spectra in a broad range of 
model parameters. 

Plane-parallel shocks with constant electron 
temperatures just discussed constitute a sub- 
class of more general plane shocks with partial 
electron heating at a shock. Electrons are also 



heated by collisions with ions downstream of the 
shock, just as in Sedov models. These shocks may 
be parametrized by the shock temperature T s , 
the postshock electron temperature T es , and the 
shock ionization age t s . If T s = T es , we arrive 
at constant-temperature shock models discussed 
above. The case with no collisionless heating at 
the shock corresponds to T es = 0, where Coulomb 
collisions gradually lead to increase in the elec- 
tron temperature T e downstream of the shock, 
which becomes equal to T s in the limit of the in- 
finite shock ionization age r s , far from the shock 
front. For partial electron heating at the shock, 
< — T es /T s < 1, just like for Sedov models. 

Multitemperature, plane-parallel shock models 
with unequal electron and ion temperatures at the 
shock should generally provide better fits to Sedov 
models. We find that this is indeed true, in par- 
ticular, we can obtain perfect fits to Sedov models 
in Figure 9. For T es = 0.5T S = 1.30 keV, the best 
plane-parallel shock fit gives T s — 5.57 keV, T es = 
0.82 keV, Ts/2 = 1.08 x 10 11 cm" 3 s, and emission 
weighted electron temperature T e = 0.37T S = 2.04 
keV. For T es = 0, we get T s = 3.36 keV, T es = 0, 
r s /2 = 9.66 x 10 10 cm -3 s, and emission weighted 
electrom temperature T e = 0.44T S = 1.46 keV. It 
may be seen that while T s and T es cannot be reli- 
ably determined by fitting shocks to Sedov models 
of Figure 9, the derived emission weighted T e and 
the SNR ionization age To compare well with pa- 
rameters of the Sedov model. 

The excellent agreement between plane-parallel 
shocks and Sedov models at high electron temper- 
atures, demonstrated in Figure 9, does not hold 
at low shock velocities. In Figure 10, we show 
X-ray spectra for an SNR with r\ = n\E = 10 50 
ergs cm~ 3 and equal ion and electron tempera- 
tures, at various times during its evolution. The 
best fit spectra produced by a plane-parallel, con- 
stant electron temperature shock are also shown. 
The fits are perfect early in the remnant's evolu- 
tion, at high shock temperatures, but they dete- 
riorate as the blast wave decelerates. Particularly 
large differences are seen at high energies, where 
the shock model spectrum is too soft. This is not 
suprising because high energy emission in Sedov 
models is produced at significantly higher temper- 
atures than at the shock, as can be seen from the 
electron temperature distribution function in Fig- 
ure 1 . Constant temperature plane parallel shocks 
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provide a poor approximation to this multitcm- 
perature distribution function, and should not be 
generally used to fit X-ray spectra of low temper- 
ature Sedov models in the 0.1 - 10 keV spectral 
range. 
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Fig. 10. — X-ray spectra for Sedov models with 
r\ = 10 50 ergs cm~ 6 and with equal ion and 
electron temperatures, for clarity displaced from 
each other by a factor of 10 in the vertical direc- 
tion. Curves are labeled by the postshock temper- 
ature (in keV). The best fit shock models (dott- 
ted curves) compare well with Sedov models (solid 
curves) at high temperatures, but they differ sig- 
nificantly at low temperatures. 



3.4. Multitemperature NEI Models 

An inadequate representation of the distribu- 
tion functions is the reason why single ioniza- 
tion timescale NEI models in general, and plane- 
parallel shocks at low shock speeds, fail to repro- 
duce X-ray spectra of Sedov models. In the latter, 
it is the presence of multitemperature plasma in 
Sedov models (see Figures 1 and 7) which is re- 
sponsible for this failure. It is possible to construct 
two-temperature, three-temperature, and gener- 
ally multi-zone "simple" NEI models which pro- 
vide a good approximation to X-ray spectra of Se- 
dov models, provided that these models are con- 
sistent with the three distribution functions dis- 
cussed previously. We constructed such models 
succesfully, and found that at most ~ 10 tem- 
perature zones are needed for Sedov models at 
low temperatures, while just a couple might be 
sufficient at intermediate and high temperatures. 
While such multitemperature NEI models might 
possibly be useful for computational reasons, they 
have been designed as a method of successive ap- 
proximation (more temperature zones give a bet- 
ter model), and not by a physical simplicity. For 
this reason we advocate the use of full Sedov mod- 
els, and not multi-zone models, in modeling of 
spatially-integrated X-ray spectra of SNRs, be- 
cause Sedov models represent a simple physical sit- 
uation of a spherical blast wave propagating into 
a uniform ambient medium. However, multi-zone 
models most likely will be very useful in spatially- 
resolved X-ray spectral modeling of SNRs. 

4. DETERMINATION OF MODEL PA- 
RAMETERS FROM X-RAY SPEC- 
TRA 

There are three parameters of Sedov models 
which may be determined from spatially inte- 
grated X-ray spectra: postshock temperature T s , 
postshock electron temperature T es , and the ion- 
ization age To (plus chemical abundances which 
we do not discuss here). It is generally difficult to 
determine all three parameters, particularly when 
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observations are available only in a limited X-ray 
spectral range or when deviations from the Sedov 
dynamics are present. We will first consider two 
limiting cases of young and hot, and much older 
and cooler SNRs. 

When young and hot remnants are observed in 
the 0.1 - 10 keV spectral range, such as provided 
by the Chandra CCDs, most of the emission comes 
from the bulk of the X-ray emitting material, with 
emission measures between 0.1 - 1.0, located be- 
hind the blast wave. Because of the low ioniza- 
tion timescales characteristic of young remnants, 
the electron temperature distribution function is 
generally flat and has a similar shape for models 
with various amounts of electron heating in this 
range of emission measures (see Figures 1 and 2). 
The ionization timescale averaged electron tem- 
perature distributions are also similar (Figure 7), 
while the ionization timescale distributions are of 
course identical and reasonably well approximated 
by a straight line (Figure 4). We have already 
shown that under these conditions plane parallel 
shocks with constant electron temperature provide 
an excellent approximation to Sedov models in the 
0.1 - 10 keV energy range. This demonstrates that 
the observed X-ray spectrum depends mostly on 
the mean electron temperature in the postshock 
region and on the ionization timescale To- It is dif- 
ficult to determine T s and T es separately without 
additional information. Observations at higher en- 
ergies, which would probe interior regions where 
the temperature distribution functions differ be- 
tween the models (see Figure 1), are not likely 
to succeed because of expected severe deviations 
from the Sedov dynamics in the interior regions of 
young remnants. 

It is possible to determine independently the 
shock speeds in young remnants through high 
spectral resolution observations. The postshock 
radial velocity of the bulk of the X-ray emitting 
material is equal to 0.75w s , which results in a ve- 
locity difference of l.5v s between the front and 
back section of the shell at the remnant's cen- 
ter. With a sufficient spectral resolution, Doppler- 
shifted spectral lines originating in these two shell 
sections could be separated, providing us with 
the missing information about the shock speed. 
This straightforward kinematic measurement is 
currently hard to accomplish because of the lack 
of X-ray instrumentation with a sufficiently high 



spectral resolution for spatially-extended sources 
such as SNRs. Such measurements would provide 
us with the shock temperature, which in combina- 
tion with the analysis described above would allow 
for a complete determination of Sedov model pa- 
rameters. Note that the high spectral resolution 
would also provide additional electron tempera- 
ture and ionization timescale diagnostics based on 
individual line ratios. 

Additional information about postshock ve- 
locities and postshock electron temperatures can 
also be independently found from optical and UV 
observations of nonradiative, Balmer-dominated 
shocks. The width of the broad Ha component 
seen in these shocks is a direct measure of the 
postshock proton temperature. In combination 
with X-ray observations, this is enough informa- 
tion to determine the shock speed in young rem- 
nants. In addition, UV and optical lines in such 
shocks are produced by electron impact excita- 
tion, (although in the fastest shocks proton ex- 
citation might also be important), and carefully 
chosen line ratios can be good diagnostics of post- 
shock electron temperatures. This is a poweful 
method for determining both temperatures from 
optical observations alone, although it is limited 
to shocks with a nonnegligible population of neu- 
tral H atoms entering the shock front. But for 
such fast nonradiative shocks one can find a com- 
plete set of shock parameters from optical, UV, 
and X-ray observations. 

In contrast to hot and young SNRs just dis- 
cussed, low shock temperatures and large ioniza- 
tion timescales in old remnants lead to efficient 
Coulomb heating and ion-electron equilibration. 
Departures from equilibration are important only 
in the immediate vicinity of the shock and in 
the remnant's interior (e.g., see top of Figure 3). 
Spatially-integrated X-ray spectra for models with 
the same postshock temperature T s but different 
T es differ significantly only at high energies, as 
shown in Figure 5 for a model with kT s = 0.49 
keV and n = 8.88 x 10 51 ergs cm~ 6 . The post- 
shock temperature T s and the ionization timescale 
To (or parameter n) can be obviously found from 
the shape of the X-ray spectrum at low ener- 
gies. But determination of the postshock electron 
temperature T es must rely on the high-energy X- 
ray spectrum. The high-energy spectrum is pro- 
duced by hot gas in the remnant's interior which 
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was shocked early in the evolution of the rem- 
nant when departures from the Sedov dynamics 
could be significant. In this situation, it might 
be difficult to distinguish between effects of colli- 
sionless electron heating and effects of departures 
from the Sedov dynamics, particularly for rem- 
nants with low shock temperatures T s and large 
ionization timescales To- For such remnants, deter- 
mination of both T s and T es is best accomplished 
through spatially-resolved X-ray, UV, and optical 
spectroscopy in the narrow region immediately be- 
hind the shock front (e. g., see Ghavamian 1999; 
Ghavamian et al. 2000). 

The postshock electron temperature T es and 
the postshock temperature T s can be simultane- 
ously determined from spatially-integrated spec- 
tra for remnants with intermediate shock tempera- 
tures and low ionization timescales To- We demon- 
strate this in Figure 11, where we show calculated 
X-ray spectra for Sedov models with the postshock 
temperature T s — 0.6 keV, with three values of T es 
(0 keV, 0.3 keV, and 0.6 keV), and for To ranging 
from 3 x 10 10 cm~ 3 s to 3 x 10 12 cur 3 s. (Again, 
the calculated spectra were folded through the 
spectral response of the back-illuminated ACIS 
CCD detector.) Note that the plane-parallel shock 
model provides a poor approximation to these 
spectra because of the relatively low (0.6 keV) 
postshock temperature. X-ray spectra for models 
with low values of To show strong dependence on 
T es , particularly at high energies, because of the 
lack of ion-electron equilibration. For such rem- 
nants, a simultaneous determination of T s and T es 
from spatially-integrated X-ray spectra should be 
possible. As the ionization timescale To increases, 
Coulomb heating becomes effective, and the effects 
of collisionless electron heating decrease in impor- 
tance. For large ionization timescales To, spectra 
in Figure 11 depend little on T es , and determina- 
tion of T es becomes difficult. 

5. SUMMARY 

We discussed spatially-integrated X-ray spectra 
of SNRs based on Sedov dynamics. Their spec- 
tral shape depends on elemental abundances and 
on three parameters: postshock temperature T s , 
postshock electron temperature T es , and ioniza- 
tion timescale t (or parameter 77). These three 
parameters can in principle be determined by com- 
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Fig. 11. — X-ray spectra for Sedov models with 
T s = 0.6 keV and with various amounts of col- 
lisionless heating: full electron-ion equipartition 
(solid curve), no collisionless heating of electrons 
at the shock (dashed curve), and partial electron 
heating at the shock with T es = 0.5T S = 0.3 keV 
(dotted curve). Models are displaced from each 
other by a factor of 10 in the vertical direction 
and they are labeled by the ionization timescale 
t (in cm~ 3 s). 
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paring models with observations for any remnant, 
but frequently only one of the postshock temper- 
atures can be found. Only the postshock electron 
temperature T es can be usually determined for 
young, hot remnants, and only the shock temper- 
ature T s can be found for old, much colder SNRs. 
In order to find both T s and T es for such SNRs, ad- 
ditional information about the shock is necessary, 
such as available from kinematic studies or from 
UV and optical observations. It should be possible 
to determine both T s and T es from the spatially- 
integrated spectrum alone for SNRs with interme- 
diate temperatures and ionization timescales. 

We succesfully explained differences between X- 
ray spectra for Sedov models with different pa- 
rameters in terms of three distribution functions: 
the electron temperature distribution, the ioniza- 
tion timescale distribution, and the ionization- 
timescale averaged electron distribution. We then 
considered a number of simple models, such as 
a single ionization timescale model and plane- 
parallel shock models, and compared their dis- 
tribution functions and X-ray spectra with Se- 
dov models. A frequently used single ionization 
timescale model was found to be an inappropriate 
approximation to Sedov models. In most situa- 
tions, plane-parallel shocks provide a much better 
approximation to Sedov models, and both plane 
shocks and Sedov models should be routinely used 
in analysis of SNR X-ray spectra. 

The NEI models described in this work, while 
constituting a basic set of simplest thermal mod- 
els, are not sufficient for interpretation of SNR 
spectra. A single plane parallel or a spherical 
shock is a crude approximation to the complex 
morphology of SNRs seen in X-ray images. For 
complex morphologies, observed X-ray spectra are 
most likely a superposition of shocks with various 
velocities, resulting in distribution functions which 
are poorly described by a single shock. In general, 
we do not expect that the NEI models presented 
here can provide perfect or even statistically ac- 
ceptable fits to SNR spectra, because of deviations 
from an idealized plane shock or a Sedov model. 
This is indeed the case for DA 530 (Landecker et 
al. 1999), 3C397 (Safi-Harb et al. 2000) and RCW 
86 (Borkowski et al. 2000), SNRs whose X-ray 
spectra were modeled with plane shocks and Se- 
dov models. More effort is clearly necessary in 
order to take full advantage of X-ray observations 



of SNRs. Improved modeling must involve a su- 
perposition of shocks with various speeds or de- 
termination of distribution functions directly from 
observed X-ray spectra, a challenging task because 
of the modest spectral resolution of CCD detectors 
onboard Chandra or XMM-Newt on. 

We based our X-ray emission calculations on 
Sedov- Taylor dynamics, in which shocks arc not 
modified by cosmic rays. But if collisionless shocks 
were strongly modified in SNRs (Jones & Ellison 
1991), then their dynamics would be poorly de- 
scribed by the Sedov- Taylor self-similar solution, 
and the resulting X-ray emission could be very dif- 
ferent than in standard Sedov models. Dorfi & 
Bohringer (1993) and more recently Decourchclle, 
Ellison, & Ballet (2000) showed that the cosmic- 
ray modification of shock structure can indeed 
have a dramatic effect on X-ray spectra of SNRs. 
This is not suprising because a significant amount 
of energy is deposited into relativistic particles 
in cosmic-ray dominated shocks, which results in 
lower temperatures of thermal gas. In several 
models calculated by Berezhko et al. (1994, 1996), 
the energy fraction contained in cosmic rays may 
be as high as 80% of the total SNR kinetic energy. 
We found that in such extreme cosmic-ray dom- 
inated SNRs the thermal gas temperature is too 
low for production of substantial amounts of ther- 
mal X-ray emission, except early in the remnant's 
evolution when its dynamics and X-ray emission 
are strongly affected by SN ejecta. Futhermore, we 
generally found fiat temperature profiles in such 
models, because of a higher fraction of energy de- 
posited into cosmic rays at higher shock veloci- 
ties. Such approximately isothermal SNRs differ 
greatly from Sedov models, which generally con- 
tain the hottest gas in their interiors. In terms 
of our distribution functions, the electron tem- 
perature distribution is flat in these isothermal 
SNRs, the ionization timescale is larger because 
of the increased postshock compression ratio, but 
the shape of the ionization timescale distribution 
is similar to that seen in Figure 4. Because of adia- 
batic expansion, the ionization-timescale averaged 
temperature is qualitatively the same as depicted 
in Figures 6 and 7. 

These trends should be also present in other 
cosmic-ray modified SNR models with lower 
particle acceleration efficiencies. A lower ef- 
ficiency is actually possible because the ideal- 



16 



ized, spherically-symmetric models calculated by 
Berezhko ct al. assume a quasiparallel shock ev- 
erywhere along its periphery, in which particle 
acceleration is most efficient. But realistic SNR 
blast waves must consist of a mixture of quasipar- 
allel and quasiperpendicular shocks, and in the 
absence of strong turbulence, quasi-perpendicular 
shocks may not accelerate particles as efficiently, 
ft is then likely that "only" ~ 50% of the SNR ki- 
netic energy is transferred to cosmic rays (Evgeny 
Berezhko, private communication), in which case 
thermal gas will produce appreciable amounts of 
X-ray emission during the late adiabatic evolu- 
tionary stage of interest here. But such SNRs 
must have significantly different distribution func- 
tions than standard Scdov models, so that they 
must produce significantly different X-ray spec- 
tra. A considerable caution is then required when 
modeling and interpreting X-ray observations of 
SNRs. The distribution functions presented by 
us are valid only for Sedov models, without any 
shock modification by cosmic rays, ft would be 
worthwhile to carry out an extensive investigation 
of the distribution functions expected in various 
cosmic-ray dominated SNR models, but this task 
clearly demands a separate effort. Note that, un- 
like for Scdov models, X-ray spectra produced 
by our plane-parallel shocks may be used even 
for cosmic-ray modified shocks. The only differ- 
ence is in interpretation of the shock temperature 
T s , whose relationship with the shock velocity 
v s is now provided not by the standard Rankinc- 
Hugoniot shock jump conditions, but by jump con- 
ditions appropriate for a collisionlcss shock with a 
particular cosmic-ray acceleration efficiency. 

Another issue connected with the cosmic ray 
acceleration in SNRs is the question of the elec- 
tron velocity distribution. In our X-ray emission 
modeling we assumed that thermal electrons can 
be well described by a Maxwellian velocity dis- 
tribution. This is most likely a good assump- 
tion for the bulk of electrons, usually referred 
to as thermal electrons, but the presence of ra- 
dio synchrotron emission in SNRs demonstrates 
that there is a high energy tail to this distri- 
bution. The most energetic (~ 10 - 100 TcV) 
electrons produce nonthermal synchrotron X-ray 
emission, as has been recently inferred for SN 1006 
(Koyama et al. 1995), Cas A (Allen et al. 1997), 
G347.5-0.2 (Slane et al. 1999), and RCW 86 



(Borkowski ct al. 2000), and modeled in detail 
by Reynolds (1998). In addition, less energetic 
electrons in the high energy tail, with energies of 
several tens of keV, can also produce lines and 
continua, through the same atomic processes as 
thermal electrons. (The continuum produced by 
such electrons is usually referred as nonthermal 
bremsstrahlung.) If these nonthermal processes 
are important in an SNR, then interpretation of 
its X-ray spectrum must rely on a combination of 
thermal and nonthermal models. A good example 
is SNR RCW 86 where a complex mixture of non- 
thermal synchrotron and multitcmperature ther- 
mal emission is present (Borkowski et al. 2000). 
It is not clear at this time how important non- 
thermal bremsstrahlung is in SNRs, because of 
the predominance of nonthermal synchrotron X- 
ray emission. But if nonthermal bremsstrahlung 
were important, energetic electrons producing it 
would also excite line emission, possibly leading 
to changes in line strength ratios and intensities. 
The possible presence of such electrons in SNRs 
is clearly an important issue for interpretation of 
SNR X-ray spectra. One might hope that with 
the high quality X-ray observations of SNRs be- 
coming available, we will be able to determine the 
distribution functions discussed in this work, and 
unambiguously separate thermal and nonthermal 
emission in SNRs. 



Reliable modeling of SNR spectra would not 
be possible without Fe L-shell atomic data, which 
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We would like to thank Andy Szymkowiak, Keith 
Arnaud, Tim Kallman, and many other people 
at Laboratory for High Energy Astrophysics at 
NASA/Goddard Space Flight Center who have 
provided continuous support for this project. We 
thank Evgeny Berezhko for discussions about 
cosmic-ray dominated SNR models and for mak- 
ing results of his calculations available to us. We 
thank Andrej Bykov for stimulating discussions 
about collisionlcss electron heating. We acknowl- 
edge numerous discussions with Don Ellison about 
cosmic-ray modified shocks. Support for this work 
was provided by NASA under grant NAG-7153, 
and by NSF through a Research Experiences for 
Undergraduates award administered by the Amer- 
ican Astronomical Society. 



17 



REFERENCES 

Allen, G. E., ct al. 1997, ApJ, 487, L97 

Arnaud, K. A. 1996, in Astronomical Data Anal- 
ysis and Systems V, eds. G.Jacoby & J.Barnes, 
ASP Conf. Series, v. 101, 17 

Arnaud, M., & Raymond, J. 1992, ApJ, 398, 394 

Beigman, I. L., Shevelko, V. P., & Tawara, H. 
1996, Physica Scripta, 53, 534 

Berezhko, E. G., Yclshin, V. K., & Ksenofontov, 
L. T. 1994, APh, 2, 215 

Berezhko, E. G., Yclshin, V. K., & Ksenofontov, 
L. T. 1996, ARep, 40, 155 

Borkowski, K. J., Sarazin, C L., & Blondin, J. M. 
1994, ApJ, 429, 710 

Borkowski, K. J., Rho, J., Reynolds, S. P., & Dyer, 
K. K. 2000, ApJ, submitted 

Bykov, A. M., & Uvarov, Yu. A. 1999, J. of Ex- 
perimental and Theoretical Physics, 88, 465 

Callaway, J. 1994, At. Data Nucl. Data Tables, 57, 
9 

Cox, D. P., & Anderson, P. R. 1982, ApJ, 253, 268 

Decourchelle, A., Ellison, D. C, & Ballet, J. 2000, 
ApJ (Letters), in press 

Dorfi, E. A., & Bohringcr, H. 1993, A&A, 273, 251 

Ghavamian, P. 1999, PH.D Thesis, Rice Univer- 
sity 

Ghavamian, P., Raymond, J., Smith, R. C, & 
Hartigan, P. 2000, ApJ, in press 

Gronenschild, E. H. B. M., Mewe, R., 1982, ApJ, 
48, 305 

Hahn, Y. 1993, J. Quant. Spcctrosc. Ra- 
diat. Transfer, 49, 81 

Hamilton, A. J. S., Sarazin, C. L., Chevalier, R. 
A. 1983, ApJS, 51,115 

Hamilton, A. J. S., & Sarazin, C. L. 1984, ApJ, 
284, 601 

Hughes, J. P., Hayashi, I., & Koyama, K. 1998, 
ApJ, 505, 732 



Hughes, J. P., & Hclfand, D. J. 1985, ApJ, 291, 
544 

Itikawa, Y., Hara, S., Kato, T., Nakazaki, S., 
Pindzola, M. S., & Crandall, D. H. 1985, 
At. Data Nucl. Data Tables, 33, 149 

Itoh, H. 1978, PASJ, 30, 489 

Itoh, H. 1979, PASJ, 31, 541 

Jones, F. C, & Ellison, D. C. 1991, Space Sci. 
Rev., 58, 259 

Kaastra, J. S., & Janscn, F. A. 1993, A&AS, 97, 
873 

Kaastra, J. S., & Mewe, R. 1993, A&AS, 97, 443 

Kato, T., & Nakazaki, S. 1989, At. Data 
Nucl. Data Tables, 42, 313 

Keenan, F. P., McCann, S. M., & Kingston, A. E. 
1987, Physica Scripta, 35, 432 

Koyama, K., Petre, R., Gotthclf, E. V., Hwang, 
U., Matsura, M., Ozaki, M., Holt, & S. S. 1995, 
Nature 378, 255 

Laming, J. M., Raymond, J. C, McLaughlin, B. 
M., & Blair, W. P. 1996, ApJ, 472, 267 

Landecker, T. L., Routledge, D., Reynolds, S. P., 
Smegal, R. J., Borkowski, K. J., & Seward, F. 
D. 1999, ApJ, 527, 866 

Lennon, M. A., Bell, K. L., Gilbody, H. B., 
Hughes, J. G., Kingston, A. E., Murray, M. J., 
& Smith, F. J. 1988, J. Phys. Chcm. Ref. Data, 
17, 1285 

Liedahl, D. A., Osterheld, A. L., & Goldstein, 
W. H. 1995, ApJ, 438, L115 

Lyerly, W. J., Reynolds, S. P., Borkowski, K. J., 
& Blondin, J. M. 1997, BAAS, 29, 839 

Masai, K. 1984, Ap&SS, 98, 367 

Masai, K. 1994, ApJ, 437, 770 

Mewe, R., Kaastra, J. S., & Liedahl, D. A. 1995, 
Legacy, 6, 16 

Pequignot, D., Petitjean, P., & Boisson, C. 1991, 
A&A, 251, 680 



18 



Reynolds, S. P. 1998, ApJ, 493, 375 



Safi-Harb, S., Petre, R., Arnaud, K. A., Kcohane, 
J. W., Borkowski, K. J., Dyer, K. K., Reynolds, 
S. P., & Hughes, J. P. 2000, ApJ, in press 

Sampson, D. H., Goett, S. J., & Clark, R. E. H. 
1983, At. Data Nucl. Data Tables, 28, 467 

Sedov, L. I. 1946, Prikl. Mat. Mekh., 10(2), 241 

Seely, J. F., Feldman, U., & Safronova, U. I. 1986, 
ApJ, 304, 838 

Slane, P. , Gaensler, B. M., Dame, T. M., Hughes, 
J. P., Plucinsky, P. P. & Green, A. 1999, ApJ, 
525, 357 

Taylor, G. I. 1950, Proc. Roy. Soc. A, 201, 159 

Teng, H., Xu, Z., & Zhang, W. 1994a, Physica 
Scripta, 49, 463 

Teng, H., Xu, Z., & Zhang, W. 1994b, Physica 
Scripta, 49, 696 



This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.0. 



19 



